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ABSTRACT 
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are analyzed to determine the storm-induced flow at various 
ocean depths, determine the associated energy increase and 
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myekrtial @MOrTLTOnN is calculated. The energy increase and 
decay is shown to vary with depth. 

An 2mbedded mixed-layer oc@an circulation model (Adamec 
e+ al, 1981) is forced with an idealized storm translating 
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I. INTRODUCTION 


Several studies of the response of the ocean to hurri- 
cane passage have been made by such authors as Fisher 
(1958), lLeipper (1967), Wright (1969), Pudov et al (1979) 
apa Fedorov et al (1979). A review of the hydrographic, 


2. se@a-surface temperature (SST), Salinity temperature 


it 


depth (STD), and expendable bathythermograph (XBT), sur- 
veys of several of these authors is given in Table I fron 
Price (1981). The majority of these observations concen- 
trate on the asymmetrical SST response, and the axpected 
ocean mixed-layer response to hurricane osassage of vertical 
mixing as the hurricane approaches, followed by upwelling in 
the wake of the storm (see e.q. Leipper, 1967; Friese, 
love eee Add tionally, oscillations on the order of the local 
znertial period have been noted in tha temperature and 
current fields in the wake of open-ocean storms. i er 
tial response becomes complicated in the shoal waters of 
continental shelves and the ability ¢t«95 detect inertial 
moticn is severely restricted (Mayer et al, 1981). The hur- 
ricane-induced upwelling decreases the mixed layer depth. 
The depth reaches a local minimum about one half inertial 


period after eye passage and continues to oscillate at about 


12 








the local inertial period for some period after the storn 
passage. Geisler (1970) presented the relationship between 
hurricane translation speed and the oscillatory response of 
the thermocline. If the hurricane translation speed is 
greater than the internal phase speed, typically about 2 m 
s-t, and if its horizontal scal2 is comparable +o the inter- 
Bat Fossby radius, currents throughout the affected area are 
controlled by a balance between the centrifugal and Coriolis 


accelerations. 


TABLE I 


Hydrographic Studies of the Sea Surface 


Temperature Response 
BOmiUemt Canes (PELGe, 1951). 











Hurricane 
Average __ central 
Un pressure A SST,,,; 








Study: Hurricane Method: Region (ms-') ~ (mb) °C) Position of * A SSTnax 
Leipper (1967): Hilda (1964) extensive post-hurricane 3 930 -§ pattern is generally 
hydrographic survey: Gulf of unclear, may be 50 km 
Mexico . to left of track (Fig. 8) 
Fedorov ef al. (1979): Extensive pre- and post-hurricane 6 980 -2 30 km to right (Fig. 3, 
Ella (1968) XBT survey: mid-Atlantic same as this Fig. la) 
Pudov er al. (1979): Tess extensive post-hurricane STD survey: 6 940 -4 75 km to right (Fig. i 
(1975) mid-Pacific same as this Fig. 2a) 
Wright (1969): Shirtey (1965) 1 pre-, | post-hurricane XBT section: 13 935 -3 20 km to night (Fig. 4) 
vicinity of the Kuroshio 
Jordan (1964): extensive pre- and post-hurncane 16 920 -2 150 km to nght 
Wanda (1956) SST reports from ships of 18 915 ~—1 50 km to night 
Clara (1955) opportunity: mid-Pacific (Figs. 2 and 3) 


ZEstimates made by Price (1981) from their figures noted 
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The response *o Hurricane Sloise as her eye vassed over 


the National Oceanic and Atmospheric Administration's Data 
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Buoy Office (NDBO) EB-10 buoy has been studied extensively 
(see €.g]. Martin, 1982; Price, 1981; Black and Withee, 
1976). EB-10 was located in the central Gulf of Mexico and 
collected the first open-ocean data under a hurricane as 
reported by Withee and Johnson (1976). Mayer et al (1981) 
reported a study of the passage of Hurricane Belle over the 
Com =amen=al shelf of the New York Bight on 10 August 1976. 
These continental shelf responses differ from the open-ocean 
respenses due to the topegraphic influences of the shelf and 
large gradients in the physical properties between the shelf 
water and the deep ocean water. 

Hurricane Frederic passed through the Gulf of Mexico 
between 11 and 13 September, 1979. The satellite data 
depicted in Figs. 1, 2, and 3 together with whatever XBT 
data collected during hurricane passage are usually the only 
type of data available to study these geovohysically, 
socially and economically important events. However, during 
+his period the Naval Oceanographic Office (NAVOCEANO), NSTL 
Station, Bay St. Louis, Mississippi soerated three data 
buoys which collected a unigue set of oC2an temperature and 
current data. This thesis studies th? ocean current 


response to the passage of Hurricane Fraieric as it passed 


14 





near the three data buoys. Additionally a comparison of the 
real data to the results proviied by the three-dimensional 
ocean model of Adamec et al (1981) is made. The hurricane 
forcing of the model is idealized rath2ar than being formu- 
lated to represent the actual hurricane, but a translation 
speed of 7.5 m s-1, the same as for Frederic, is used. 

The buoy current data providei by NAVOCEANO indicate 
that the flow associated with the passage of Hurricanes 
Frederic hada large inertial component after the storn 
passed. The rate at which the inertial flow damps, and the 
rate at which it propagates with depth, are determined from 


the buoy data and compared to the model. 
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A. BACK GROUND 

Ten current meter array records collected at three sites 
on the continental shelf of the southeastern United States 
(Fig. 4) during the passage of Hurricane Prederic were pro- 
vided by NAVOCEANO. Newdepecteod intrag. 5, the arrays coa- 
Sisted of Aanderaa spar buoys anchored +) the bottom. The 
anchor cable held instruments at three depths at each cf the 
Stations one and two and at four depths at station three. 
The shallowest instrument depth was 19 m at station two, 
while the deepest was 457 m at station three. The bathyne- 
try (Fig. 5) shows all three stations ara on the continental 
shelf. Data were provided for a centinuous period from a 
few days before the storm passage, about 2100 Greenwich Mean 
Time (GMT) 12 September 1979, until a few days after «he 
StOrIM passage. 

faeoredeptn at Station one is only about 100 neters. 
Ppeetcemmoe?) reported difficulties with Meters in a similar 
coastal environment because of the presence of very strong 


horizontal gradients and strong topegraphic effects on the 
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AI Deep Site 
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Menet Ie 


Hegre &. gest position track 9£ Hurricanes Frederic and 
Ga2e2Ons Of CUEDTent Meter arfays (provided by 
pie oNeay, NAVOCEANO) . 


currents. Mayer et al reported intense, first-mode, near- 
Prete period oscillations in water of gabout 70 meter 
depth. Similar observations at sites in shallower water of 
about 50 meter depth resulted in only weak, heavily damped, 
second-node oscillations. 

The current meter records used in this thesis also show 


€ Variation in response among the stations with different 








bot+com depths. At station three where the water depth is 
about 465 meters, an increase in current nagnitude starting 
about 1920 GMT 12 September 1979 is very apparent in the raw 
data record (22g. Ga). This recori also shows an 
oscillation which cycles through one wavelength about once 
every 24 hours (1440 minutes). Pius o eae decays 
smoothly over about a week, and at first glance it appears 
Sommave a period of about 22 hours. A more detailed analy- 
Sis cf the record shows the zaro down-crossings to be very 
nearly 1450 minutes apart during ‘the largest oscillation in 


the record. The average period for sevan cycles during this 


interval is 1430 minutes (23.3 hours) between down-cross- 
ings. The local inertial period at station three is 1473 
minutes. Shortly after the surface response, the rapid 


increase in current magnitude is also obvious at the depth 
of 457 meters at station three (see Fig. 6b). The average 
pericd of the oscillation at this greater depth, 1370 min- 
utes or 22.8 hours, is less than the period néar the sur- 
face, but the damping of the oscillation is not obvious at 


457 meters. At station one where the oottom depth is only 


fai 


apout 100 meters, ‘the increas? in current speed is detect 


ble in the raw data record near «he surface (Fig. 78), but 
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any period of oscillation onthe order of 24 hours or a 
damping rate is difficult to detect. An increase in current 
velocity on 12 September 1979 near the bottom at station one 
(Fig. 765) is not apparent and it can be seen that current 
magnitudes are larger on 9 September than they are on 12 
September. | 

To determine if the observed currents were extraordi- 
nary, some idea of the likely currents in the area of obser- 
vations is helpful. A dominant feature of the Gulf of 
Mexico surface currents is the Loop Current (Fig. 3). This 
Msea Clockwise current of about 50 to 2039 cm s~*t and 90 ¢o 
150 km width (Leipper, 1970) which entsrs the Gulf in the 
west as the Yucatan Current, and exits through the Florida 
Straits as part of the Gulf Stream system. The Yucatan 
Current flows north from Honduras between the Yucatan 


Peninsula and Cuba into the central east Gulf and forms the 


western section of the Loop. FUSener noOwrn the current 
flows east, southeast, and ¢ast again. Elsewhere in the 
northeastern Gulf, the surfaces flow is generally cyclonic 


along the coast until turning southward near Louisiana. 
maecomeare knewn to have detached Erom the Loop Curyen< 


(Leipper, 1970) and could move ints the area of «he 
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instrument arrays. Leipper (1970) shows a systematic devel- 
opment and breakdown of the Loop Current which places the 
MeLtnein edge of the current further north in the spring 
than in other seasons (Fig. 9). Molinari (1978) suggests 
the northward intrusion of the Loop Current is not seasonal, 
but that climatological results are biased by temporal san- 
pling techniques. Hurlburt and Thompson (1980) presented a 
numerical study of Loop Current intrusions and eddy shed- 
ang. They obtained theoretical expressions for the eddy 
diameter and penetration distance of the Loop Currant into 
the Gulf and concluded that if only vorticity dynamics are 
considered, tne interaction between the Florida Shelf topog- 
raphy and the pressure field results in a balance which 
stops the northward penetration of the Loop Current. All of 
the current array stations are further north than 29 deg 
North, and although the waters of the Loop rarely extend 
Breese fas north, Huh st al (1981) examined an intrusion of 
Loop waters as far north as the data buoys which collected 


+he raw data used in this thesis. Lewis Hoe poss up] sO 


ry) 


Pommamoguorl whether or not the Loop currcent is in the area 


of the current meter arrays during the peso of 


observations, but it seems unlikely. 
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Figure 8. Eyp2 cal Surface currents in the Gulf of Mexico 
Or the months of oe August, September. The 
eneral location of the ee arcayS 1s indicated 
y the box. (U.S. Navy, 1965). 


Another factor which might be expectsd is the influence 
of tcpography. The principle cf conservation of potential 
vorticity causes the current t> follow the bathymetric con- 
Sours. The proximity cf the buoys to the coast suggests 
that near-shore property gradients between the shelf water 
and the deep water may have important dynamic and thermnody- 
namic consequences. The DeScts Canyon is also in the area 
cf cbservations aidmoeomey ¢cend to enhance cross-shelt 
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Figure 9. The spring intrusion of 1966 as indicated by 
overlays of the 150-meter contour lines from the 


BUNCE pee of all spring cruises of 1966. 
May 1964 and November and Decenber 1965 are 
1970) 0 for supplemantal information (Leloper, 


B. RAW DATA 

The NAVOCEANO data consisted of north-south (v) and 
sast-west (u) velocity components and temperature (T) az 10 
minute intervals. Temperatures above 21.5 9C ware not 


recorded due to thermistor limitations. Temperatures above 


Zo 








21.5 °C occurred only near the surface and complete tempera- 
ture records were available only from the deeper 
instruments. Also, NAVOCEANOD reported apparent internal 
clock synchronization problems in the current meter record 
at 64 meters depth, station one. 

The hurricane passed the data buoys with a translation 
speed of about 7.5 m s~!, and as depicted in Fig. 4, the 2ye 
passed about 100 km to the west of station three at 2200 GMT 
12 September 1979 (Julian day 255). Maximum winds of 115 
MmeoseCGGClGsed at a radius of about 30 km from the storm cen- 
ter. The inertial periods for the data buoy locations range 
fromy 25 hours, 57 minutes at station on2 to 24 hours, 33 
Minutes at station three. That is, the inertial period at 
all of the stations is approximately one day. Sincs? the 
diurnal +tidal period is 24.8 hours ani the semi-diurnal 
mada lem=ticdsas 12.4 hours, di#ficulties could arise in sep- 
arating any inertial motion from tidal motion before the 
storm arrival. A spectral analysis would most likely noi 
have sufficient resoluticn to separate inertial and diurnal 
Won iG The first harmonic of the inertial frequency would 
also be indistingquishable from the semi~iiurnal ‘«idal fre- 


quency. However, the large increase in currents, and «hus 
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kinetic energy, observable in the raw data records would be 
easily detected in an energy spectrun. 

The initial data records provided by NAVOCEANO span 
non-coincident periods. However, all rscords contained data 
from 1920 GMT 07 September 1979 to 1820 GMT 21 September 
ero. The raw data records (Figs. 6 and 7) show an obvious 
respcense to ere huseacane forcing starting about 12 
September 1979. This hurricane forcing results in increased 
currents, with maximum values of 130 cm s7~! near the sur- 
baCe.. Comparison of Pigs. 6a and 6b shows that the forcing 
response was transmitted as deep as 457 neters in less than 
cne half an inertial period. The response to hurricane pas- 
sage is more difficult to detect in most of the raw data 
records of stations one and two and it is therefore not pos- 
sible to estimate the rate at which the response is verti- 
cally propagated with depth merely by looking at the raw 
data reccrd. The ability to detect the increase of energy 
associated with storm passage at any station and at any 
level is discussed in the section ‘titled "Procedures for 
determining relative energies". Although this thesis does 
not analyze the temperature records directly, +the response 


+o hurricane passage is also apparent in the temperature 
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mecord (Fig. 10). The temperature record shows an 
escillation in the thermocline and an apparent cooling at 
the 437 m depth at station three, followed by a gradual 


Warming to the end of the period. 


DEG C 


STATION 3 DEPTH 937 MEERS TEMP 





Seed 2 LO bi l2 13) «14 > POSueet 7 weGouerld. (cGl Se) 
DATES DENOTE 1920 GMT 


Figure 10. Raw temperature record for 7-21 September 1979 
at a depth of 437 m at station 3. 


Oe LATA MANIPULATION 

The initial records were shortened to include only the 
pericd 1920 GMT 7 September 1979 to 1820 GMT 18 Sepnptember 
oo. This period was the longest possible period included 
by all records and was a logical choice since the hurricans 


eye passed all three buoys between 2100 GMT 12 September and 
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0300 GMT 13 September 1979. By examining plots of the data 
(eid. Figs. 6 and 7) it was determinei that significant 
storm effects with some pre-storm and som2 post-storm data 
would be included in the shortened record. 

ProgresSive vector diagrams (PYD's) were plotted for 
each record during the above period. These diagrams indi- 
cated that the inertial motion was superposed ona mean 
current which was quite different at each location. The 
mean current was not of primary interest and attempts were 
made to remove the mean current from the records so that the 
magnitude, damping period, and associated energy of the 
inertial motion could be determined. In the first attempt, 
the mean current over the first three inertial periods 
(pre-storm passage) was subtracted from all of «he u veloc- 
ity ceomponen+ts andthe v velocity components. The three 
inertial period average started at 1920 GMT 07 September, 
which is five days before the storm center passed the near- 
est buoy. Comparison of the initial PYD"s (Fig. 11) and 
Similar diagrams of the same data with the three inertiai 
period average removed (Fig. 12) was maie. 

THemenictal FyO'S indicate 3 general -Eransport in direc- 


sions which, in most cases, are nearly the same before and 
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after the storm-induced perturbation. The perturbations are 
genezaliy reminiscent of inertial motion, i.e. oscillations 
With a period equal to the local inertial period. PVD's of 
purely inertial motion would transcribe a simple circle for 
each inertial period. The circles would overlap if more 
than one inertial period were plotted. Where the first 
three days of a record were similar to th2 mean flow, remov- 
ing the thre2 inertial period average from the record 
resulted in more nearly overlaoping circular PVD's. Renov- 
ing the three inertial period average was not effective in 
isolating the inertial component of the flow when the first 
three days of the record were not obviously in the general 
transport direction. Therefore, the mean current was chang- 
ing significantly during/following the storm at many o= the 
current meter locations. A more effective method of isolat- 
ing the inertial component from each of the records was 


required. 


D. ISOLATING THE INERTIAL MOTION 

In a further attempt *o isslate the large inertial con- 
Penencmecna. waS ObViowts in the original records, 4 funning 
average cver the inertial period was computed each three 


home = ond =nen subtracted from the initial PVD recoris. The 








averaging technique is necessarily applied at the central 
point of each inertial period, and as a result the half- 
periods at the beginning and end of the records were lost. 
However, the transient response is not great at these times 
and the information lost is inconsequential. Next, the 
average inertial period record is linearly interpolated 
between 3-h values to obtain th same nunber of data points 
as the initial PVD data. PVD's of th2 running inertial 
period average data (e.g. Fig. 13) show that the majority of 
the inertial flow in the raw data is removed by the running 
average. The running inertial period averaged PVD's show 
the general transport direction as well as some features 
which should be noticeable in 3 model of the ocean response. 
A surge to the northeast, onshore, is noticeable in thea 
meccrd of Fig. 13a, and there is Some indication that the 
post-storm current is greater than pre-storm, i.e. the hori- 
zontal displacement between inertial period marks is greater 
aiter the storm surge than before. The best example of 
+hese features in any of the records 1s shown in Fig. 13b. 
Actiegemcou-ne weet, LOollowed by a counterclockwise isop back 
into the pre-storm flow direction is evidsnt. Additionally, 


the post-storm flow is of much greater nagnitude than tine 
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pre-storm flow. These diagrams were obviously a better rep- 
resentation of the mean flow throughout the record than the 
previously computed three inertial period averages. These 
PVD's also show that the mean current was indeed influenced 
by the bottom topography. Comparison of Fig. 4 and repre- 
sentative FYD's at stations one, two and three (Figs. 14 and 
15) show that the mean flow follows tha bathymetric con- 
cours. The current moves rapidly toward the northwest ax 
station one during and for some time following the period of 
hurricane passage (Pig. 14a). A counterclockwise loop which 
would be associated with inertial motion is not seen as it 
Moran the staticn two record {(Fig. 1%b}. Comparison of 
Figs. 13b and 15 shows the flow is along th2 contours but in 
opposites directions at the 251 and the 437 or 457 meter 
depths. Also, a shoreward surge similar to that at the 
near-surface level of station three (Fig. 13a) is seen at 
least one inertial period later near the bottom (Fig. 15). 
item cendency §2O return to the pre-stomm flow direction is 
much slower at «he near-bottom records and is not completes 
by the end of the record. 

The running inertial-period averages are then subtracted 


rom the initial PYD data point-by-point and PYD's of the 
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Figure 13. 
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Bosulieant Gaza are plotted (Fig. 16). Examination of these 
resultant PVD's shows that station thr2e exhibits the most 
obvious inertial response, possibly due to less coastal or 
bottom influence. The resultant PVD's for stations one and 
cwio 6€6(nOot shown) also show large inertial responses. The 
station one record still exhibits a strong flow which is 
non-inertial. The station two recori shows that the 
majority of the non-inertial motion is removed. 

It 1s not possible from these plots to determine vre- 
cisely when the response to the hurricane forcing first 
begins at each depth. Since the rate at which the energy is 
transmitted vertically is of interest, thea next procedure is 
to calculate the energy associated with each inertial period 


at all available levels. 


Ee PROCEDURES FOR DETERMINING RELATIVE ENERGIES 

The PVD's which were computed in the above steps indi- 
cate the average of the current record is not zero. Also, a 
breakdown of these records into inertial periods indicate 
this nonzero mean current which does exist is not constant 
eaeeougnou-. the record. An interactive computer program was 
developed to remeve the mean current from a record consist- 


ing of inertial rlow superposed on a mean current. The nean 
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velocity is calculated as the distance between the starting 
and ending points of the record divided by the duration of 
the record. An equal portion of this mean current is sub- 
tracted from each ten minute interval of the record. The 
mean current is removed from the record one inertial period 
at a time rather than for an entire record since the mean 
current 1s not constant throughout the entire record. 
feemrecle Sor Tadaus R was Fit to each of the inertial 
pericd records. The circle coordinates (x,y) are computed 


using a mean velocity. Given that 


x R cos (9) 


y = 8 sin(O) 
The mean velocity is computed as the time derivative of the 
position as: 

u = dx/dt = -R sin(@) dOsit 

v = dy/dt = R cos(m) de@sit 
where 

d@sdt = -27/(N at) 

and N represents the number of points us2d to make the cir- 
ene. Successive x and y coordinates of the circle are then 
computed as: 


Coe =x ete (uy it ) 


y1 yO + (v dct) 
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me Starting  spoSition is x = 0. anduye= — | owe cee Vee 270 
degrees of a 360 degree circle centered at (0.,0.), and the 
successive positions are created in a tlockwise fashion. 
The radius of the circle is adjusted using an error mini- 
mization program. This program calculates the distances 
between the positions required to create the circle and the 
positions used to compute the PVD's, and selects the circle 
which gives the smallest sum of positional distances. One 
fourth of the maximum dimension of the plotted data is used 
as the initial radius for the minimization subroutine. 
While the radius is held constant, the sum of the difference 
between an “average error" and an “individual position 
error" is computed and divided by the number of differences. 
This error difference is used in the following calculations. 
The distances between all points of the computed circle and 
the data are first compared to locate the data point and the 
Circle point which are closest. These points are used as 
the starting positions for the distance error computations. 
The “average error" is the average of the distances between 
*he consecutive points of the jata and the computed circle. 
The “individual position error" is the distance between a 


point on the circle and a data point. The error difference 
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for cne radius is stored and, the radius is slowly incre- 
mented to twice the original magnitude, storing the error 
difference at each increment. The radius corresponding to 
the minimum value of error diffarence is used to create the 
circle which best fits the plotted data. 

The circle created using this adjusted radius is used as 
the best fit. It is then possible to estimate the energy 
per unit mass associated with each inertial period. The 
inertial velocity (V=f£R) associated with the circle radius & 
and coriolis parameter ££ is computed. Onewnalreot “the 
square of this velccity is assumed to be the energy per unit 
mass. The citcle fitting error difference computed as above 
is divided by the radius of the fitted circle and referred 
Semcommeamcc mmr ventS error. The energy vaiues are plotted 
versus the inertial period at each leval at each station 
(See Figs. 17 and 18). The radius for sath inertial period, 
the relative RMS error of fitting the circle and the mean 
current removed from each inertial period are listed in 
Table II for the records from station one at a depth of 49 o 
and from station three at a depth of 21 on. Also, cheem’s.- 
Dezathoene program is tested using exact circles as input 


data. The resultant errors for given radii are also listed 
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in Table II. The exact circle error shows the error is near 
zero when the motion is circular. The 2rror is found to be 
@ smaller percentage of the radius, i.e. the relative RMS 
error is smaller, at station three than at station one. 
This supports the previous results of the motion being more 
nearly circular, and thus more inertial, at station three 
than at station one. 

A good fit between the positions of the generated circle 
and the positions from the real data with the inertial mean 
removed is obtained for the records of station three where 
the water is deepest and the coastal influence is least 


(Fig. 19a). The agreement in positions between the two 


curves at station one is generally less since these records 


are affected more by the botton, and by tidal motion which 
is not necessarily removed. Obviously non-inertial motion 
remained in the station one record (Fig. 19b). Example 


plots of the record after all means were removed are shown 
inet 1 dias 20 « The figures show nearly overlapping circular 
moticn except for different radii. Plots of the energy val- 
ues versus inertial periods are given in Fig. 17 for the 
level nearest the surface at 2ach of the three stations. 


Peecayeevalues of 4000 at station 3 correspond to an 
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true circle to the manipulated data accounts for some of 
this difference. When the data points ar2 not circular, the 
circle is adjusted +o some value less than the largest 
dimension of the plotted data. When the data form a circle, 
*he largest dimension of the plotted data and the dimension 
Gieetne Circle nearly coincide. Por stations two and three, 
Set Peaincrease in kinetic energy is indicated in Fig. 17 
apenge the fifth inertial pariod. These energy values 
remain high for two to three inertial periods before rapidly 
decreasing toward the pre-storn levels. Sit eas Lots EOE 
the next deepest level at 2ach station are shown in 
Pgs 183. t stations one and two the energy peaks rapidly, 
falls off rapidly and then increases to the end of the 
mecord. + station three, a slower increise is noticed, and 
is fcllowed by a nearly constant value to the end of the 
record. A slower increase in anergy is found (Fig. 18b) at 
the 437 mand 457 n levels at station three. The energy 
values do net appear +o have peaked by the end of the 


EaCOrd. Samce "the =tecords of Fig. “18 are within 20 m of 
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(iD 


each other, it would be r2asonabdl 
would be nearly identical. The difference could be due to 


bottom reflection of energy since the 457 m instrument was 
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positioned only 8 meters from the bottom. Awe the very 
least, the procedures used to determine the energy associ- 
ated with the inertial motion should produce results accu- 
GBaee to within a factor of two and should also indicate the 
trends. No confidence is plac2i in the details of the fluc- 
tuating values. However, the large pzaks and long tern 
increases shown in the energy records are thought to be 


Srgn2fscant. 


9000. 


Se . 


tL PUM 


Av) 
cS 
) 
ec 


(By 
CS 
i) 





Nie Clse Clete Oss 


Figure 17. Energy (m2 s~*) versus ine 
neaz surface levels at sta 
three. 


49 








dud. 








sU0. 179 m Sta 2 
al 
ue. 
Be 64 m Sta l 
Lal m ota 3 
mt) Z 
| 
le. ——e 
Q en U 5 
Me ie Wels 
BOo!, 
437 m Sta 3 
B00. 
=a 
iz.) 
> yoo. — 
be 
\ 
\ 457 m Sta 3 
200. y~ Ae 
Gee I | I 
0 e UY 6. 8 10 l2 
U3/ AND U5S7 M 
Figure 18. Enera 2 s72) yersus i ms a) El edaat 
: the ag esis ot ach Station abd (b) See ero 
bottom levels at station thres. 





Figure 19. 


KILOMETERS 


KILOMETERS 


eee are real 


and 


5 


GONE INERTIAL PERIOD 
NORTH 


KILOMETERS 


ONE INERTIAL PERIOD 
NORTH 


KILOMETERS 


data record 
e for (a) inertial peri 
(b) ae eaa p2riod a’ 


31 


VW 














ARROWS INOICATE INERTIAL PERIGOS 
NORTH 





0 
Ce) 
qe 
uJ 
wi -10 
x 
© 
| 
<= 
~20 
~30 EAST 
-10 -5 O S 10 15 20 2s 
KILOMETERS 
ARROWS INDICATE INERTIAL PERIGODS 
g NORTH 
6 
U 
W 
< 
Lad 
2 
x 
SQ 
= 
= 
0 
-2 
-4 
-6 -4 -2 0 2 y 5 
KILOMETERS 
Figure 20. PVYD's after removal 9% mean currents at a depth 
2 of (a) 21m and (b) 257 m at station three. 


2 








III. THE EMBEDDED MIXED LAYER--OCEAN 
CIRCULATION NODEL 
he BACKGROUND 
Adamec et al (1981) tested an embedded mixed laver- 
oc@an circulation model. This model embedded Garwood's 


(1977) model for predicting mixed layer depth and jumps in 
*emperature and velocity at the base of the mixed layer ina 


six-level, primitive equation »pen-ocean model developed by 


Haney (1981). These authors remarked in their conclusion 
that additional tests, “including comparisons with observa- 
tional data", were being planned. This thesis constitutes 


an important step in those plans becaus2 of the availability 
of suitable observational data. The model used in the 1981 
study is expanded to simulate the threa-dimensional ocean 
Bespense to a translating hurricane-type forcing field. 
This model forcing represents an axisymmetric hurricane 
eaeansiaeang at a cOnStant speed of 7.5 a s-~?. The model 
ocean basin is a 960 km square with six layers of varying 
depth between the surface and a free-slip bottom at 400 nn. 
The horizontal resolution is 15 km on a 65 by 65 grid. The 


level depths and thicknesses are represanted in Fig. 21. 


The time step used in the 1981 study, 450 s, is increased to 


aS 





600 s. A brief explanation of the model as used in this 


thesis 1s given below. 


Be MODEL FORMULATIONS AND BOUNDARY CONDITIONS 

The governing equations ar> the primitive Navier-Stokes 
equations, the continuity equation, the hydrostatic equation 
and the equation of state (sea Adamec st al, 1981). The 
hydrestatic assumption is appiied, and the ocean is assumed 
incompressible with density being a linear function of temp- 
erature only. The coriolis force vari2s appropriate to the 
latitudes of the domain from 25°9N to 33.6°N. The change in 
coriolis parameter (f) with latitude is computed as a finite 
difference at each north-south grid point. 


There are no fluxes of mass, momentum or heat across tha 


bottem (flat) or side (vertical) boundaries. Also, the 
rigid lid approximation (w=0 at 22:0) is nade, which 
requires the vertically averaged motion be zero. Applying 


this approximation to the Navier-Stokes equations provides 
prediction equations for the vertical shear currents (see 
Adamec 2+ al, 1981). VeEvica tlyesvereseging che continuity 
equaticn over the mixed layer and applying w=0at 22:20 


yields 2 prognostic equation for the vertical velocity at 


*he base of the well-mixed layer, w(-h). The entrainment 
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Figure 21. Representation of the six mixei-layer--ocean 
circulation model layers. Computation isvels 
acum vediidetlare ji apornts agted an each level. 
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velocity is predicted by the mixed layer model. Together 
With w(-h), this ylelds a prognostic equation for mixed 
layer depth. 

The axisymmetrical portions of the hurricane forcing are 
shown in Fig. 22. The tangential and radial stresses are of 
the same form. They are calculated from a prescribed wind, 
which increases linearly from the storn 2ye wall to the 
radius of maximum winds, and then decr2ases as r-!/72 to f= 
360 km and then linearly to r = 450 km. The radius of mnaxi- 
mum winds is 45 km and the inner boundary of the eye wall of 
the storm is 4.5 kn. This wind profile results in a wind 
stress curl which is zero inside the eys wall, increasing 
linearly from zero at the eye wall toa maximum at 45 kn, 
zero from 45 to 360 km, and negative from 360 kn to the 
boundary of the wind stress (about 450 km). The maxinum 
tangential stress corresponds to a maximum wind speed of 50 
m s~t, The small value of radial stress (-12.9 4iPa max) is 
due to cross isobaric flow of about 209. 

PFelow the mixed layer, vertical diffusion is appiied to 
the momentum or temperature equations with a vertical eddy 
viscosity or eddy conductivity coefficent (both equal to 3.5 


cm s7~?2 at all depths below the surface nixed layer). In the 
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Normalized Forcing 


Radius [km] 


Figure 22. Forcing functions for the ocean model, each 
normalized by the value of the function at the 
radius of maximum winds. The tangential and 
radial stress components have the same forn 
(solid line) with maximum values of 35.9 and 
-12.9 eer respectively. The Sea se £Lex 

a mn 


acess S a maximum value of -84 

Adamec et al, 1981). 

mixed layer, however, the above formulation is not 
aeprepriate. The depth of tha mixed layer is the boundary 


between intense turbulence and the musth less turbulent 


waters beneath. 5 


C. ENTRAINMENT AND MIXED LAYER MODEL FORMULATION 

In the mixed layer, prognostic equations for the nixed- 
layer average (bulk) values cf the vertical component of 
turbulent kinetic energy andthe total turbulent kinetic 
energy are derived using the bulk, sactond-order closure 
methcds of Garwood (1977). Computation of the entratznament 


buoyancy flux allows calculation of the jJownward fluxes of 
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heat and momentum associated with entrainment at «he base of 
the mixed layer. 

The above mixed-layer formulation allows the fraction of 
turbulent kinetic energy available for mee nd at the base of 
+he mixed layer to be other than a constant. This turbulent 
kinetic energy is diagnostically dependent upon che surface 
buoyancy flux and the surface friction velocity (known 
boundary conditions). It is thus possible to couple «he 
diagnostic mixed-layer formulation with =he prognostic ocean 
circulation model in a numerically feasible fashion. 

Foth increases and decreases in mixed layer depth must 
be considered in any mixed layer model. The easier of these 
two events te formulate is the increasing mixed layer depth 
case. The method requires determination of the entrainnent 
heat flux and then imposing this heat flux on the given 
temperature profile. Added vazrtical resolution near «he 
base of the mixed layer is given sinc? the basé is not 
required to coincide with any of the prescribed model lev- 
els. This increased vertical resolution is very important 
Since the thermocline profile determines the potential 


energy of the upper ocean. The dynamics of the mixed layer 


‘2 


and the ocean circulation are dependent o this potential 


energy. 
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Formation of a new, Shallower mixed layer occurs 
whenever the turbulent mixing is not able to penetrate all 
the way to the previously established depth of the mixed 
layer. This layer reformation occurs when warming of the 
surface layer occurs with no increase in the wind shear to 
sustain the mixing. Numerical formulation of this event is 
difficult because the previous structure at the interface 
can not readily be preserved. Thus energy budget problems 
may arise when deepening to the prior mixed layer depth does 
occur. This model uses a numerical procedure (see Adamec et 
eu, 6©6: 198 1) that preserves potential 2nergy to ensure the 
deepening rate is as correct as possible when the layer 
again deepens +o the earlier mixed layer depth. This fea- 
“ure of the model is not tested here because the solar flux 


is zero throughout the simulation. 


D. Diner C STABILITY CONDITION 

It is assumed that the mixed layer is dynamically unsta- 
ble and the underlying water column is normally dynamically 
stable, However, dynamic instability of the underlying 
water column carn sometimes occur. In this model, vertical 
Fluxes of heat and momentum between levels are imposed so 


that the gradient Richardson number remains greater than or 


aid 








equ tO a Gritical value. This generalization of the more 
common convective adjustment is referred to by Adamec et al 


(1981) as "dynamical adjustment". 


ape OF THE DYNAMICAL AND MIXING PROCESSES IN THE 
t 


The mixed layer and the dynamic portion of the model, 
which is a level model that oredicts the average of a quanr- 
tity in a layer, are coupled in two phases. Horst, -aavec- 
tive and diffusive changes in the upper ocean are calculated 
in the dynamic part and put ints a form useable by the mixed 
layer model. Then, the surface flux and antrainment changes 
are calculated by the mixed layer model and transmitted to 
the dynamic part of the model. A special treatment cf the 
level which contains the base of the mixed layer is required 


(see Adamec et al, 1981). 


Fie COMPARISON OF THE MODEL AND THE DATA 
The model provides a four-dimensional space-time view of 
the ocean response to a moving hurricane. The data provided 


by NAVOCEANO is *«wo-dimensional, time and depth, at three 


th 


Geeeeecreoe Stations. All three locations are to the no 


rh 
p- 
Q 


jie wsroame, track beyond the tadius of maxinum winds of +h 


ty 


storm. The differences in instrument depths and the limited 
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mumber of stations do not allow the extension of the data 
into a third dimension. The comparison o£ the data and «he 
model, then, is an attempt to: 


e compare the general en) eee characteristics of the 
model results to see i hey are reasonable in light of 
the observed data results, 


e choose a_ data set from the model which is similar to 
the raw data in number, depth and Location relative to 
the storn, 


e compare this data to the NAVOCEANO data to show whether 
or not the model 1S providing a realistic response to 
the simulated hurricane forcing, and 


e examine the four-dimensional model results for charac- 
teristic ocean response to hurricane forcing which 
could not be observed at the data stations. 


G. GENERAL CHARACTERISTICS OF THE MODEL 
iicwenatraiwoesztion of the hurricane 2s at grid point 
(23, 16). Since the grid spacing is 15 km, grid point 
(23,16) corresponds to a distance in the x direction of 330¢ 
Kipeand ian the y direction of 225 kn. Model output includes 
north-south and east-west velocity components, Vea 
temperature, T ; and mixed layer depth, H. The initial val- 
ues cf these gqguanties are given in the following table. 
The guantities u, v, and T ara extracted at both three hour 
Bicermvels fOs rhe Cntire grid at all six levels and at ten 
Minute intervals for selected grid points and levels. The 


Mixed layer depth (MLD) is stored at three-hour intervals. 
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TABLE Lert 


Initial values of model output variables 


T at level 
uo H i 2 3 s 2 eC 
(om "s-1) (m) (°C) 
0 0 2.06 50). 0 29.2 23.0 Zone 18.4 lesen 


The results predicted by the moiel at 36 and 43 h at level 1 


are 


presented in Figs. 23a--d and Figs. 24a--d, respec- 


tively. Using the storm translation speei of 7.5 m s-!, the 


distance of the storm from the bottom of the grid is given 


in Table IV. The storm center has tracked off the grid by 


hour 


the 


Boe 


and 


24 and the southern boundary of the storm has passed 
edg2 by hour 48. 


TABLE IV 


icane distance from bottom 9f grid at six-hour intervals 


DISTANCES HOUR 
330 0 
492 6 
654 2 
816 18 
978 24 

1140 30 
1302 36 
1464 42 


-_ <== eee <2 cm a Se em ee ee ow oe ew aw SD ee ee Se Oe ee ee ee es ee a ee eee ee ee ee a eee eee ae eee eee ee 


Several important features are noticeable in Pigs. 23 
Zu; 


Doemecimcrent ctelds at 36 and 48 hours are "out of 
phase", i.e. the areas of, westerly currents at hour 36 
are generally easterly at hour 48. 
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nt motion is circular over n 


ta early +h2 entire 
th an apparent rotation period of about 24 


Tange from neg- 


yer wt 
positive values of 


e Current velocities in th 
ative values of about 7 


about 80 cm s~!., 


2 mixed la 
0 “@cm-t to 


e The mixed layer deepens acne teas ft TOmene Cight of 
the storm, and ranges from minimum values of about 30 oa 
to maximum values of almost 60 om. 


e A ridge of cold, upwelled water lies in the wake of the 
Storm,  Tnas ridge is shown best in the level 3 temper- 
dtire sacontour of ee Ds Heo OSCOmnoOn ons. cold 
ridge 1S an upwelling/downwelling wake with an inertial 
period, GUpwelling is_in the areas of minimum MLD's 
2 «dp gn Pig. ee and downwelling is in the area of 
maximum MLD's (58.90 in Fig. Sea Minimum MLD's are 
reached at the end of the upwelling cycle and maximun 
MLD's are reached at the and of the jlownwelling cycle. 


These results are consistent with the results derived fron 
the NAVOCEANO data, which shows that the primary ocean 
current response to hurricane passage when no mean current 
exists is circular motion at the local inertial period. The 
reversal in flow direction every 12 hours is particularly 
evident in the u velocity component records for level six 
Siown in Fag. 26. The model predictions describe the ccean 
response over the entire grid, rather than jus= at some 
selected points. Trajectories are plotted at several loca- 
tions perovendicular to the storm track (see Fig. 27). The 
Surface-layer trajectory that began on the left side of the 
track (Figs. 27a) clearly shows the surge +o the ieft as 


Siemec-OUMlapproaecnes, Strong insrtial motion following storn 


65 








passage and a net transport to the south. PVE 2a JSCtory 
along the storm path (Pig. 27b) shows a surge to the left, 
recovery toward the right and a small net displacement. The 
trajectories that began on the right (Figs. 27¢ and 27d) 
show a net deflection northward and to the right. These 
trajectories show that different transport profiles are 
expected at different locations relative to the storn. This 
feature is not obvious in the NAVOCEAND data due to the 


limited number of locations of the observations. 


H. SELECTED POSITION AND DEPTH DATA 

A data set similar to the NAVOCEANO data set was derived 
by extracting model variables at ten minute intervals at the 
three positions and depths shown in Table V. These time 
series of model variables (Figs. 28 and 29) are clotted sin- 
ilarly to the raw data plots of Figs. 6 and 7. The observed 
(Figs. 6) and predicted (Fig. 28) u-conponents show very 
Similar variations aithough the model data are much 
Smoother. Also, although the magnitude of the velocity con- 
ponent at level 2 of Fig. 28a is somewhat smaller than the 
21 m magnitude shown in Pig. 53, the magnitude ar level 6 
(300 n) of Fiq. 28b is much smaller than che 437 n magni- 


SUdeCMsown in Fig. 6b. Ti2wewplacetion is that not enough 
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Figure 25. Model temperature prediction for level thres at 
(ay Hour 36 ard (bj) hour 48. 
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energy is being vertically transmitted by the model. Bion 
+he model motion and the observed flow are inertial and nei- 
ther has completely damped by the end of nine days. The 
model u components of the velocity at levels one and six are 
exactly 180 degrees out of phase (Fig. 30a). That is, the 
velocities or the surface and bottom model layers are in the 
opposite directions. The NAVOCEANO velocities are also 
exactly out of phase immediately after storm passage, 
although they become nore in phas@® latar in the record (see 
meg. 30D). For any given level of the nodel, the u and v 


Variations are 99 degrees cut of phase. 


TABLE V 
Positions and depths for model ten minute interval data and 
corresponding NAVOCEANO data. 
DISTANCE 
INERTIAL F ROM MID NAVOCEANO 
GRID PERIOD SLORY CENTER LAYER INSTRUMENT 
POSITION MODEL NAVOCEANO MODEL NAVOCEANO DEPTH DEPTH 
meme scese sf iM 
807,33 24.50 24.54 105 100 874 49. 
75.0 64, 
150.0 92. 
235 24. 30 24.34 90 85 19.0 19, 
150.0 17 9e 
300.0 324, 
Seo 25699 230 120 Zs io a0 2 
1 ee) Zale 
15,0220 437. 
30020 457. 


Three-hour running averages by inertial period were con- 


puted for the extracted model data using the same technique 
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as used for the NAVOCEANO results. A direct comparison of 
Figs. 13, 14 and 15 to the extracted moiel results running 
inertial period averages is mad2 in Figs. 31, 32, 33. The 
model near-surface residual flows are very Similar to each 
other (Figs. 31a and 32a), as are the lower level results 
Megs. 31D, 32b and 33). Net transport to the north is 
indicated in all records of th2 model. The NAVOCEANO data 
indicate net transports that are closely tied to the local 
topography, which is absent in the model. The NAVOCEANO 
results also show opposite current directions for near-sur- 
face and mid depth. These observations provide interesting 
examples of results not shown in the model, and which 
require addtional investigation to explain. 

In order to make a more quantitative comparison between 


+he model and the NAVOCEANO inertial oscillations, the model 


(D 


results from the selected stations shown in Table V wer 
treated in the same manner as the NAVOCEANO raw data. Cir- 
cle radii, relative RMS érror and mean current removed are 
Peo camenaple Vil) similarly to Table II. Relative RMS 
errors are similar to those for station «hree (Table ITI) 
although maximum circle radii are not as large. Resultant 


calculations of energy (cm? s-@) versus inertial period are 
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plotted (Figs. 34 and 35) similarly to the energy plots of 
Figs. 17 and 18. These figures show that the major differ- 
ences in the model and the NAVOCEANO data, at least ina 
qualitative sense are in the rate of propagation of energy 
with depth. It appears that less energy is propagated ver- 
tically in the model compared to the observations. The 
energies in Figs. 35a and 35b are much smaller in comparison 
*o the values in Pig. 34 than are the similar results of the 
observed data. Also, there is no slow inzsrease in the hori- 
zontal kinetic energy near the bottom in the model as 
observed in the two deepest NAVOCEANO records. The level 
four records in Fig. 34a suggest that energy is being con- 
tinually propagated into the layer because the energy levels 
do net appreciably decay. Detailed calculations aré neces- 
sary to establish the energy flux mechanisms from the model 
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IV. SUMMARY, CONCLUSIONS AND RECOMMENDATION 


A. SUMMARY 

The NAVOCEANO data provided several physically 
significant results. The SUErents measured by the 
instruments were affected by topography and by ‘the storm. 
The flow tended to be along tha bathymetry contours except 
that the storm imparted a surg= to the flow consistant with 
wind stress associated with the storn. The station one 
records do not depict easily iientifiabl= inertial notion. 
Station one is in shallow water of about 100 meters, and 
shoreward motion to mid-depth is evident prior to the storm 
passage. An offshore flow is evident fron about the tine of 
the storm passage for a duration of about two inertial peri- 
ods. The tendency at each station is for the post-storn 
flow to return to the same direction as the pre-storm flow. 
Energy input by the storm has generally dissipated by a fac- 
tor of e71 after three or four inertial periods, and the 
post-storm flow magnitude is greater than the pre-storm mag- 
nitude. The horizontal kinetic energy assoclated with iner- 
tial motion in the near-surface water increased more than 
two orders of magnitude as the storm passed and then decayed 


over six or seven inertial periods. At the deepest 
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instrument levels, the kinetic energy associated with iner- 
tial motion increased slowly and had not peaked seven iner- 
tial periods after storm passage. The mii-level instruments 
Tecorded either a sharp peak after storm passage or aless 
rapid increase and little decay ts the end of the record 
depending upon lecation. The near-surface and néear-botton 
flows are about 180° out of phase for at least five or six 
inertial periods after storm passage. 

The mixed layer-ocean circulation model results are very 
Similar in many respects to the NAVOCEANO observations. The 
model does not include bathymetry, and therefore bathymetric 
effects are not reproduced. The model does show, however, a 
storm-induced surge, inertial motion danping over several 
inertial periods at the surface and bottom layers and dif- 
ferent magnitudes in pre-storm and post-storm flow. The 
model predictions indicate an almost instantaneous near-in- 
ertial energy propagation raté with depth. The flow at the 
surface and at depth have the same form, increasing in mnag- 
nitude and decaying at very nearly the same rate at each 
level. A+ mid levels of the model, the kinetic energy level 
remains nearly constant. The @cuwerent oscillations “in ne 


mixed layer and in the bottom layer are 180% out of phase. 
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Trajectories of the model results show that the net dis- 
placement at some stations is nearly zero while at other 
stations the net displacement is tens of kilometers. Also, 
the direction of surge and of the residual flow is dependent 
on lccation relative to the storm track. The model shows a 
zone Of upwellina/downwelling in the wake of the storm which 


has a near-inertial frequency. 


ae CONC LUSIONS 

The bathymetry of the ccean is very important in deternm- 
ing the direction of flow of currents at all levels. hes 
PemsceoOr the flow is not included in this ocean circulation 
model and must be subtracted from the real ocean data if 
consistent comparisons of model results and current observa- 
tions are to be made. Dive Clit lor mmom cho Storm 2S tomdrive 
the average motion in a direction consistent with the wind 
field of the stora. Changes in direction occur over a very 
Short time period and are observed at all depths. As the 
storm passes, inertial motion results due +0 a balance 
between the cent+rifugal and coriolis accelerations. The 
kinetic energy associated with this motion remains nearly 
constant near the surface for ‘two to three inertial periods 


and decays over an additional three to four inertial periods 
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tO a post-storm level greater than the pre-storm level. The 
available data do not show the energy levels returning <+o 
the pre-storm level before the end of the record. There is 
an okvious time lag of the storm which increases with depth. 
The exact time of the initial effect at any particular depth 
is difficult to measure since the magnitude of the velccity 
does not instantaneously increase, but increases slowly with 
the onset of the storn. Nonetheless, the downward propaga- 
tion of this energy is nct instantaneous but requires on the 
order of one-half inertial period to propagate one-half kil- 
ometer. The energy propagates to a depth of about 450 m in 
one-half of an inertial period. At this depth the energy 
increases for at least seven to eight inertial periods. The 
mixed layer-ocean circulation model closely simulates the 
inertial motion and damping rates near the surface. The 
vertical propagation of energy is not, however, very Similar 
+o the NAVOCEANO observations. The model does provide 
insight into <the three-dimensioanal ocean and is impressive 
inedaseirgu2shbing between cesults at locations varying with 


zespect to the storm center. 
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C. RECOMMENDATIONS 


The primary differences between the NAVOCEANO data and 


the model results are in the vertical propagation of kinetic 


(n 


energy associated with storm passage. Further research <i 


Oo 


required to isolate the mechanism of energy propagation s 


that it can be included in the model formulation. 
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